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Abstract. 

In recent years several Dwarf Novae (DN) systems have been ob¬ 
served in quiescence, when the accretion rate is low and the WD pho¬ 
tosphere can be directly detected. The WDs are observed to cool after 
the DN outburst from high effective temperatures to lower effective tem¬ 
peratures (Teff) thought to be indicative of the thermal state of the deep 
interior of the WD. Sion has argued that the most likely energy source 
for this quiescent luminosity is the gravitational compression of the WD 
interior, which rejuvenates an otherwise cold WD into a much hotter 
state. We are undertaking a theoretical study of the compressional heat¬ 
ing of WD’s, extending down to the very low time averaged accretion 
rates, (M) ~ yr“^, applicable to the post-turnaround CV’s 

(the “TOADS”). Nuclear burning is unstable at these (M)’s, so we have 
incorporated the recurrent heating and cooling of the WD throughout 
the classical novae limit cycle. In addition to self-consistently finding the 
range of Tgff as a function of (M) during the cycle, we also self-consistently 
find the ignition masses. Comparing these theoretical masses to the ob¬ 
served ejected masses will tell us whether the WD mass in CV’s is secu¬ 
larly increasing or decreasing. We close by comparing our results to the 
accumulated observations of quiescent DN and making predictions for the 
colors of low (M) CV’s in quiescence that are applicable to searches for 
faint CVs in the field and galactic globular clusters. 


1. Introduction 

Dwarf Novae (DN) systems contain a white dwarf (WD) accreting matter at 
time-averaged rates (M) < 10 “®Mq yr“^ from a low-mass (< O.SMq typically) 
stellar companion. At these (M)’s, the accretion disk is subject to a thermal 
instability which causes it to rapidly transfer matter onto the WD (at M {M)) 
for a few days to a week once every month to year. The orbital periods of these 
binaries are usually less than 2 hours (below the period gap), but there are also 
DN above the period gap, > 3 hours (see Shatter 1992). The M onto the WD 
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can be low enough between outbursts that the optical/UV emission is dominated 
by the internal luminosity of the WD, not the accretion disk. Recent HST/STIS 
spectroscopy has spectrally resolved the WD’s contribution to the quiescent light 
and found effective temperatures Teg ~ 10,000 — 40,000 K (see our Figure § 
and Sion 1999). 

The measured internal WD luminosity is larger than expected from an 
isolated WD of similar age (ks Gyr), indicating that it has been heated by 
accretion (Sion 1985). Compressional heating (i.e. internal gravitational energy 
release) appears to be the main driver for this re-heating (Sion 1995). Sion’s 
(1995) early estimate for internal gravitational energy release within the WD 
(of mass M and radius R) was L ~ 0.15GM{M)/R. However, we show in 
§2 that most energy is released in the accreted outer envelope, giving L ~ 
3kTc{M)/finip, where /r ~ 0.6 is the mean molecular weight of the accreted 
material, nip is the baryon mass, and k is Boltzmann’s constant. The theoretical 
challenge that we address in §3 is how to calculate the WD core temperature, 
Tc, as a function of (M), and thus find Teg. Because of the unstable nuclear 
burning and resulting classical novae cycle, the envelope mass changes with time. 
This allows the core to cool at low accumulated masses and be heated prior to 
unstable ignition. We use nova ignition to determine the maximum mass of the 
overlying freshly accreted shell, and find the steady-state (i.e. cooling equals 
heating throughout the classical novae cycle) deep interior temperature of the 
accreting WD, Tc, as a function of (M) and WD mass. 

In §4, we compare our theoretical work to STIS observations and infer (M) 
on the timescale of 10® years, critical to constraining CV evolutionary models. 
Figure Q shows that DN above the period gap are hotter than those below the 
gap, and have (M)’s consistent with that expected from traditional CV evolution 
(e.g. Howell et al. 2001), even those that involve some “hibernation” (Shara 
et al. 1986; Kolb et al. 2001). The result is more surprising if the much 
weaker magnetic braking laws of Andronov et. al. (2001) are correct. We also 
predict the minimum light {My) of (M) < yr“^ CVs in quiescence, 

allowing for discovery of the predicted large population of CVs with very low 
mass companions (< O.IMq) that are near the period minimum (Howell et al. 
1997). Observations already show that the WD fixes the quiescent colors of these 
CVs and our calculations are useful for surveys in the field that were discussed 
at this meeting (e.g. 2DF, SDSS, see Marsh et al. 2001 and Szkody et al. 2001 
contributions here), as well as HST CV searches in globular clusters. 

2. The Basic Physics of Compressional Heating 

It is important to make clear what is meant by internal gravitational energy 
release or compressional heating within the white dwarf. Most are familiar with 
the energy released {GMnip/R) when a baryon falls from a large distance to 
the stellar surface. This energy is deposited at, or near, the photosphere and is 
rapidly radiated away. This energy does not get taken into the star, as in the 
upper atmosphere (where T -C Tc) the time it takes the fluid to move inward 
is always much longer (by a factor of at least Tc/T) than the time it takes for 
heat to escape. Thus, once accretion has shut off, or diminished (such as in DN 
quiescence), this energy release is no longer relevant. What is relevant is energy 
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release deep in the WD due to compression by the freshly accreted material and 
hydrostatic nuclear burning. That energy takes a long time to exit, will still be 
visible when accretion has halted, and sets Teg in quiescence. 

Let’s start with the simplest view of compressional heating. In the non¬ 
degenerate outer atmosphere, the fluid is slowly falling down in the WD grav¬ 
itational field, g = GM/B?. A fluid element falls a distance of order the scale 
height, h = kT/grUpg, in the time it takes to replace it by accretion, giving 
L ^ {M)gh {M)kT/gulp. This exhibits the scaling that appears in the for¬ 
mal viewpoint, which is to consider the heat equation 



f)c 

= T^ +Tv-Vs 
ot 


dL 


( 1 ) 


where e^v is the nuclear burning rate, s is the entropy, and v = —{M)f/A'Kr‘^p 
is the slow downward advection speed from accretion. In pressure units, and 
neglecting nuclear burning and the time-dependent piece, this becomes 

The entropy decreases inward (i.e. the envelope and core are not convective), 
so this is an outward going L. The entropy profile is fixed by the temperature 
gradient needed to carry the luminosity outward and thus we simultaneously 
solve equation ( 2 ) (putting back in e^) with the heat transport equation, using 
opacities and conductivities from Iglesias k. Rogers (1996) and Itoh et al. (1983). 
For an analytic understanding, we integrate through the non-degenerate enve¬ 
lope, where s = klniT^^"^/p)/gnip. Here, L is nearly constant, giving T®’® oc . 
Integrating down to the isothermal core yields L ~ 3kTc{M)/ gmp. 

Now, let’s turn to the degenerate carbon/oxygen core. For the (M)’s and 
typical M = O. 6 M 0 WD of interest here, all of the entropy is in the liquid ions 
at Tc ~ 10 '^ K. The time it takes to transport heat through the interior is 
r\j lO^yr <C M/(M), so the core is isothermal and any compression is far from 
adiabatic.Q Due to uncertainty from the classical novae cycle, we don’t know 
whether the C/0 core is secularly increasing in mass, but if it were, almost all 
of the work of compression goes into increasing the electron Fermi energy. The 
integrated heat release would only he L ^ 15kTc{M) /giirip (Nomoto 1982) for 
a 0.6Mq C/0 WD, where pi ~ 14 is the ion mean molecular weight. Because of 
the mean molecular weight contrast, this is about a factor of five smaller than 
that released in the envelope. Thus, the entropy drop through the accreted layer 
is larger than that across the core, making the accreted layer the main source of 
compressional heating. 


^This is in contrast to the rapid accretion rates (M) S> 1O“®M0 yr“^ considered for more 
massive Type la progenitors, where the interior undergoes nearly adiabatic compression (see 
Bravo et al. 1996). 
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Figure 1. The Hydrogen/Helium envelope and outer core in tem¬ 
perature and pressure for three different values of accumulated mass, 

Mace = 0.5 X 1O“^M0, 1.5 X 10"^Mo and 2.5 x IQ-^M© for M = O.GMq 
and (M) = lO“^^M0/yr. The external surface luminosity of the WD 
and the corresponding is also listed. The part of the star off to the 
right of the figure {P > 10^^ dyne cm“^) is the isothermal inner core. 

3. Finding the Equilibrinm Core Temperature 

For our initial study, we dropped the time-dependent term in equation ([I|), and 
presumed that the C/0 WD mass was constant throughout the classical nova 
cycle, thus only accounting for the compressional heating and ctv in the H/He 
layer. This method improves on that of Iben et al. (1992) by allowing the 
accreted envelope mass to change through the 10® year classical nova cycle. 
Early in the cycle, the mass of the accreted layer is small, compressional heating 
is small, and the WD cools. Later in the cycle, the accreted layer becomes thick 
enough that compressional heating along with slow hydrogen burning releases 
a sufficient amount of energy to heat the core. As the WD has a large heat 
capacity, reaching the equilibrium Tc where the heat exchanged between the 
envelope and core averages to zero over a single classical nova cycle takes ~ 10® 
years. Since this time is usually shorter than the time over which (M) changes, 
we construct such equilibrium accretors for a given M and (M). 

To do this, we first fix Tc at the outer edge of the C/0 core at a pressure 
high enough so that the changing accumulated mass has little direct effect. With 
a radiative outer boundary condition, we then integrate our structure equations 
with equation (|^) to find the envelope state for an (M) and accreted layer mass. 
See Figure for examples of the resulting T-P relations. We then evaluate the 
luminosity across the chosen location (the right edge of the plot in Figure 2) 
for a number of different accreted layer masses up to the unstable ignition. The 
ignition mass is found by comparing the T and p at the base of the accreted 
(hydrogen rich) layer with the analytic ignition curves in Fujimoto (1982). 
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Figure 2. Luminosity at the outer edge of the core as a function of 
accumulated mass up to the classical nova ignition for the equilibrium 
model of Figure Positive is outgoing and epochs of core cooling and 
heating are indicated. 


We vary Tc to find an equilibrium model, where the “core luminosity” (Lcore) 
averages to zero over the classical nova cycle as shown in Figure The quiescent 
Tefj which is expected for the same cycle is shown in Figure ^ At the nova out¬ 
burst we assume that the accreted shell is expelled, and that due to the rapidity 
of this event, it does not appreciably heat the WD. The resulting equilibrium 
core temperatures for (M) = 10“^°MQ/yr are Tc = 9 x 10® K, 7.5 x 10® K and 
8.5 X 10® K for M = 0.4,0.6 and I.OMq. The M = 0.4Mq star is hotter than the 
M = 0.6Mq star because it has a larger maximum accumulated mass that leads 
to a longer period of core heating. For a O.GMq WD, the core temperatures are 
rc/10®K = 4, 5.3,12.2 and 18.0 for (M)/M© yr'^ = 10"^^, 3.2x10“^^, 4.2x10“^® 
and 10“®. This also gives us the Teff range during the classical nova cycle, which 
we now compare to observations. 


4. Comparison with Observations 

In contrast to Sion’s (1995) simulations which only directly addressed the imme¬ 
diate cooling after outburst, our work has focused on the release of heat from the 
deep interior. The long-timescale nature of heating in these deep layers allows a 
comparison between our results and the observed Teff’s. For CVs in the field, the 
large set of STIS observations by Szkody et al. (2001) and previous observations 
(Urban et al. 2000) provide detailed spectra of quiescent WDs in DN. Figure 
^ shows the current Tgfj measurements and compares them to our results. The 
measurements are made during deep quiescence when the accretion luminosity 
is negligible and are intended to be long enough after the outbursts that other 
emission mechanisms (e.g. Pringle’s (1988) suggestion of radiative illumination 
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Figure 3. Effective temperature of the WD as a function of accu¬ 
mulated mass for equilibrium model of Figure |^. 
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Figure 4. The quiescent WD Teg’s in DN as a function of binary 
orbital period. The open circles are Cycle 7/8 STIS observations from 
Szkody et. al. (priv. commun.), while the solid circles are from the 
compilation of Urban et al. (2000). The single point for GW Lib is 
from Szkody et al. (2000), which is cold enough (Teg ^ 11,000 K; Van 
Zyl 1998) to be a ZZ Ceti pulsator in quiescence. The thick vertical 
lines are our current predictions for the Tgg range for the stated masses 
and (M)’s. The thin solid line is the ZZ Ceti instability strip (Bergeron 
et al. 1995). 
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of the WD) have faded. In that case, we are observing heat directly from the 
deep WD interior. The number of such systems will increase due to upcoming 
all-sky surveys (such as SDSS and 2DF, see Marsh et ah 2001), and will push 
to lower (M) systems with long quiescent intervals. 

We can already give some general insight from the initial results shown. 
Our work indicates that below the period gap, (M) ~ Mq yr“^ and the 
WD masses are in the range O.G-I.OMq. This agrees with the expectation from 
Kolb &: Baraffe (1999), who find (M) Ri 5 x yr“^ at an orbital period 

of 2 hours presuming angular momentum losses from gravitational waves alone. 
Above the period gap, the Tefj is higher, and we estimate (M) ~ IO^^Mq yr“^. 

We predict that a O.GMq WD above the gap has a core temperature of 
1.8 X 10^ K and, if in equilibrium below the gap, Tc = 7.5 x 10®K. An interesting 
question is whether the WD will have time to cool this much as it traverses 
the gap. If not, then the WD’s below the period gap will be hotter than our 
calculation implies. We estimate this cooling time from the current WD cooling 
law (e.g. Chabrier et ah 2000) along with the heat capacity of the core. 


Acooi ~ 10 


0 


1.8 X lO^K 


2.5 


A 

dt 


SksT 

mrup 


M 


(3) 


which gives At ~ 0.5 Gyr. Since this is comparable to the estimated time spent 
in the gap (Howell et al. 2001), our equilibrium assumption below the gap is 
likely safe. However, note that about 0.2 Gyrs after accretion halts, the WD 
will enter the ZZ-Ceti instability strip! 

These results also have bearing on the search for faint CVs in globular 
clusters. The expected population might well contain many low (M) systems 
that spend much of their time in quiescence. These CVs are commonly searched 
for via the presence of hydrogen emission lines or X-ray emission (as recent 
Chandra observations have found; Grindlay et al. 2001a, Grindlay et al. 2001b), 
and this method is fruitful. However, for the very-low (M) systems, the disks 
can be very dim and the quiescent X-ray emission too faint even for Chandra. 
We now show that these systems (as well as CVs crossing the period gap or those 
“hibernating” post-novae, Shara et al. 1986) can be identified by their position 
in a color-magnitude diagram (CMD). By using our theory of the thermal state 
of the WD, it is possible to predict the broadband colors of CV systems without 
dependence on the disk luminosity. 

An excellent example is NGC 6397 (King et al. 1998; Taylor et al. 2001). 
Figure shows a CMD of NGC 6397 with our initial results that provide a 
relationship between Teg and (M). The lines were produced by superposing a 
WD with the maximum Tgff for the indicated (M) with a MS star. Except for 
near the WD cooling line (dashed curve), where the WD becomes completely 
dominant, the I magnitude is set by the MS companion. The large dots along 
the 10 “®Mq yr“^ and 10 “^®Mq yr“^ lines indicate where the MS companion is 
0.3, 0.2, 0.15 and 0.1 M©, and two additional points at 0.09 and 0.085 Mq are 
indicated on the W~^^Mq yr“^ line. This immediately provides a number of 
candidate systems (namely, those data residing in this part of the CMD). 

The circled points are the “non-flickerers” (Cool et. al. 1998) recently 
reported by Taylor et al. (2001). The three at I ~ 22.25 are very strong Ha 
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Figure 5. Color-Magnitude Diagram of NGC 6397. The dots are 
members of the cluster that are below the MS (HST observations by 
King et al. 1998) and the O’s are the “Non-Flickerers” from Taylor et 
al. (2001). The MS line is from Baraffe et al. (1997) for [M/H]= —2.0 
and an age of 10 Gyr (actual cluster [M/H]= —1.95). The dashed line 
is from Bergeron et al. (1995) for DA WDs with log 5 = 8. The lines 
connecting the MS to the WD sequence are our current calculations of 
the WD-I-MS binary at the specified (M). The highest T^g during the 
classical nova cycle has been used for the WD in each case (see Figure 
No disk has been included, which is a safe guess for the low (M) 
systems. All curves have been put at the distance and reddening of the 
cluster, (m — M)j = 12.05 and E{V — I) = 0.288. 


absorbers (consistent with a DA WD) and were not detected by Chandra (Grind- 
lay et al. 2001b). These authors had discussed these systems as possible helium 
WDs with millisecond pulsar companions, though, given our work, we would 
claim that these are hot WDs with ~ O.15M0 MS companions. In addition, 
the population of data points in this diagram with respect to our theoretical 
curves will eventually constrain CV evolutionary scenarios. If we assume many 
of the data points are CVs, we already see that most systems with high (M) 
have 0.15 — 0.3Mq companions. The large number of data points below the 
(M) = 10 “^^Mq yr“^ line could well be the long-sought post-turnaround sys¬ 
tems with (M) = 1O“^^M0 yr“^ and companion masses < O.O9M0 (Howell et 
al. 1997). 
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5. Conclusion and Future Directions 

We have evaluated the action of compressional heating of accreting WD interiors. 
Most of the compressional energy release takes place in the accreted envelope, 
and is thermally communicated to the core. The maximum envelope mass is set 
by the unstable nuclear burning that causes a classical nova runaway and most 
likely expels the accreted mass. We have constructed equilibrium accretors which 
have constant core temperatures such that the heat lost from the core when the 
envelope is thin (i.e. right after the classical nova) is balanced by that regained 
when the envelope is thick. This equilibrium determines the Teff of the WD 
throughout the classical nova cycle. Our models agree with the observations of 
Dwarf Novae in deep quiescence and imply (M) ~ Mq yr“^ just below the 
period gap and (M) IO^^Mq yr“^ just above the period gap for WD masses 
in the range O.G-I.OMq. 

Though our initial efforts have met with apparent success, there is still 
much to be done. First we need to investigate the relevant parameter ranges 
(e.g. M and (M)) within our initial scheme. The most critical parameter to 
vary is the metallicity of the accreted material, lowering to values relevant for 
globular cluster science. We must also survey our initial assumptions: 

• WD excavation or accretion. The assumption that there was no net mass 
loss or gained by the WD through the classical nova cycle must be relaxed. 
This will allow for cooling of the WD due to adiabatic expansion after mass 
loss, or heating if it is increasing in mass. The large C/0 fractions seen 
in some novae ejecta (Gehrz et al. 1998) might indicate that the WD is 
decreasing in mass at (M). 

• A self consistent accounting for the ignition masses, including varying the 
metallicity. Our initial work used Fujimoto’s (1982) results for simplicity, 
but these are limited at low (M)’s. A careful comparison to more modern 
nova calculations (e.g. Prialnik & Kovetz 1995) must be carried out. 

• Thermal evolution of the WD. We also need to concern ourselves with the 
secular change of (M) due to the decreasing companion mass or changing 
angular momentum losses (such as a drop in magnetic braking at the 
period gap). When this occurs on a timescale comparable to the WD 
thermal time, it is possible that the WD will not reach the steady-state 
solution we have assumed. This memory of previous higher (M) epochs 
could well allow the WD to be hotter than the equilibrium accretor. 

While these steps are unlikely to change our understanding of the com¬ 
pressional heating mechanism, each is essential for applying our work to the 
observations. As assumptions are relaxed and investigated, much more will be 
learned about both the state of CV systems and their evolution. 
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